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Nonperturbative effects in event shape distributions can be characterized by shape functions derived in the eikonal 
approximation or, equivalently, from soft-collinear effective theory. The use of energy flow operators and the boost 
invariance of the Wilson hnes of soft gluons in the shape functions leads to a proof of universality for power corrections 
to the mean values of event shapes, without invoking the single gluon approximation. 

In this Comment, we elaborate the formahsm of soft-cohinear effective theory (SCET) presented in for the shape 
functions for thrust, C-parameter, and angularities. We show that several well-known universality relations for mean 



, ^ values of these observables can be derived without invoking the approximation of single soft gluon emission. These 



relations can be derived either in the language of SCET at leading order in the expansion parameter A = -^Aqcd/Q 
0,01 or, equivalently, in the eikonal approximation 



The central result described in Q i^ an operator exp7;ssion for the nonperturbative shape function for the distri- 
bution in the event shape e (where e — > in the two-jet limit): 

"S- 5e(e) = -^ = ^TrV|(X„|r„Ffi|0)|'5(e-e(X„)), (1) 
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where do is the total cross-section at leading-order in as, Nc is the number of colors, the sum is over final states 



Xu in the ultrasoft sector, and YmYn are Wilson lines of ultrasoft gluons corresponding the outgoing quark and 

antiquark in the n, n light-cone directions. The delta function fixes the contribution to the event shape e{Xu) for 

'■^^ . each final state in the ultrasoft sector equal to a value e. We will develop below a way to eliminate the explicit 
M-i. 

I , sum over these states. 

The event shapes we will consider are the C-parameter and angularities (which include the thrust). It will be 
• useful to re-express these event shapes in terms of the pseudorapidities ry^ of the final state particles, where 



77,=lncot-^, (2) 
■ measured with respect to the thrust axis t. For massless particles, these are equivalent to the rapidities 

with the z-axis taken to be along t. In terms of the pseudorapidities, the angularities 0] and C-parameter fj] can 
be expressed as: 



Q^t^coshry, 

In this form, the relations between the event shapes implied by invariance under Lorentz boosts will be made more 
manifest. 

A striking prediction from the analysis of these event shapes |^-|l2j is the universality in power corrections to 
their mean values, 

(e) = (e)pT + Cey^, (6) 
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with A a universal parameter and Ce a calculable coefficient that depends on the observable, 

2 

cc = Stt. 



1-a' 



(7) 



(We limit our attention to angularities for a < 1.) In an approximation that we discuss below, their effect also is to 
shift the event shape distributions. 



d<j 



PT 



— > -5-(e - CeA) 
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These relations were derived in Refs. 0,0,0] from the assumption of a "dispersive" representation for as{pt?) 
considered as an analytic function of the scale /i. In Refs. 14, thev were abstracted directly from the form of 
resummed perturbation theory. A more general approach in [5|, |y, llll | replaces the shift by a convolution with the 
shape function defined as above. In 0| and j^l the role of energy flow was explored in a manner closely related to 
our discussion below. 

While we have seen how the single gluon approximation may be implemented in the effective theory to reproduce 
the well-known relations between the shifts to event shape distributions, the shape function already contains within 
itself all the information needed to prove Eq. © for the mean values, without making any further approximations. 
We simply need a few tools to unveil this information. We will concentrate here on the shape function associated 
with ultrasoft radiation. Its factorization from perturbative collinear radiation is discussed in the context of SCET 
in Q IE 1^ 1 or equivalently in jsj . 

Consider first the angularity distribution. In this case, the shape function is 



^Tr^(0|ytyJ|Xj(X.|r„y^|0)Wr,-i ^ |k,^|e-l'"l(i-'^)j. 



(9) 



We now define a transverse energy fiow operator, by analogy to the energy flow operator introduced in Ref. [l^. 
defined through its action on a state N: 



(10) 



where the sum is over the particles i in state N . Making use of this operator, we can write the shape function as 

(11) 



SrAra) = ^Tr^ (0|y,Ii;t<5 {^-^^j rf^e-l''l(i-'')fT(r?)) (XJ y„F« |0) . 



The operators in the matrix element no longer contain any reference to the final state X„, so we may perform the 
sum over intermediate ultrasoft states, leaving 



SrAr,) = ^ Tr (0| yIy^5 {^a - ^ j d^c-^^^'^'-'^'^ St{^)^ Y^Y ^ |0) , 



(12) 



Now consider the properties of this matrix element under a Lorentz boost along the z-direction with a rapidity t]' . 
Inserting factors of U{k{r())'^U{k{r()) = 1, 

SrM = ^ Tr (0| U^UyIy,Iu^US (^Ta - ^ J drje-\^\^'-'^£T{v)^ U^UY^J^U^U |0) . (13) 
The vacuum |0) is invariant under Lorentz boosts. Wilson lines are also invariant: 



r„(0) = Pexp 



^9 



dsn ■ Aus{ns) 
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where a — e ^ , as n an and n a ^n. (This is also known in SCET as type-Ill reparametrization invariance 
The only change is in the operator £t{v)'- 



= UiA{r,')) J2 Hl^iv i) \N(A-'h)) 



ieN 



Y^\ki\5{Tl + i -r^,)\N{h)) 



(15) 



ST{v + ri')\N) 



Thus, we conclude that 



= J_Tr(0|Fkt^ {ra-^j dr;e-l''l(i-'^)£:T(r; + r;')) r„F^ |0) . 



(16) 



This result does not imply that the leading power correction simply shifts the argument of the perturbative event 
shape distributions, as the delta function (which is just the shape function to leading order in perturbation theory) is a 
highly nonlinear function of the energy flow operator and sits sandwiched between Wilson lines in the matrix element. 
If we neglect correlations between these operators, however, we derive a delta function for the shape function, and 
reproduce the shift in the distribution, Eq. (0) 5, 12]. 

However, consider the implication of this formula for the first moment of the Tq distribution. Taylor expanding-'^ 
the delta function in Eq. (|16|l , 

SrAra)^5{Ta)~5\Ta)^ jdr^e-\^\^^-'^^^TT{Q\YU^£T{7i + v')^^^^ (17) 

From this relation we see that, remarkably, the first moment is completely independent of the value of the rapidity 
appearing in the energy flow operator £'7^(7/), as we are free to choose any value for 77'. Thus, we may take the matrix 
element containing the £t operator out of the integral over r/, leaving us to evaluate simply: 

dr^Q-H(i-a) ^ (18) 
3 1 — a 

For the C-parameter distribution, the whole argument would be the same, except that the rapidity integral appearing 
here would be 

■X 

drj — = 37r. (19) 

3 cosh 77 

Thus, the shape functions satisfy 

SrAra) = S{Ta)-S'{Ta)-^A+--- (20) 

1 — a 

SciC) = 6iC)-S'iC)ST:A+-- - , (21) 

and the universal quantity A is given by 

-4 = Tr (0| YlY,l£T{n)YnY^ |0) , (22) 

which is in fact independent of the value of 77 chosen in £x{ri). When convoluted with the perturbative distribution, 
Sr^{Ta) and SciC) reproduce the universality relations of Eq. © for the first moments of the distributions. We have 
thus established these results without appealing to the one-gluon approximation. 



^This expansion is valid in a region of size A near the endpoint, where Aqqq <g; A <g; Q, so that smearing the distribution over this 
region yields an expansion in powers of Aqcd/A- 
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